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According to the unified model, the differences observed in different classes of AGN are due 
to anisotropic obscuration. This implies that in the hard X-rays, at energies where the radia- 
tion can pierce the obscuring material, the average spectral characteristics of different classes of 
AGN should be the same. Here we present a study of the 17-250 keV INTEGRAL IBIS/ISGRI 
average spectra of Seyfert galaxies, aimed to constrain similarities and possible differences be- 
tween different classes of AGN. The sample we consists of all the 172 Seyfert galaxies detected 
by INTEGRAL IBIS/ISGRI at z<0.2 during its first 7 years of observations. Of these, 45 are 
Seyfert Is, 31 Seyfert 1.5s, 67 Compton-thin Seyfert 2s, 1 1 Narrow Line Seyfert Is, 1 1 Compton- 
thick Seyfert 2s, and 9 LINERs. Comparing the average hard X-ray spectra of the different classes, 
we found that, while the spectra of Seyfert Is, Seyfert 1.5 and NLSls are consistent, Compton- 
thin Seyfert 2 show a clear excess in the 20-60 keV band over the spectrum of Seyfert Is and 
Seyfert 1.5s. This excess might be due to an intrinsically stronger reflection component or to a 
clumpy structure of the absorber in Compton-thin Seyfert 2s. 
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1. Introduction 

Active Galactic Nuclei (AGN) are thought to be powered by accretion onto a super mas- 
sive black hole (Rees 1984). Seyfert galaxies host AGN, and are commonly divided in Seyfert Is 
and Seyfert2s. While Seyfert 1 galaxies show both broad (FWHM > l,000kms _1 ) and nar- 
row (FWHM ~ 300— l,000kms _1 ) optical emission lines, Seyfert 2s show only narrow lines. 
Seyfert 1.5s show an optical spectrum intermediate between those of Seyfert Is and Seyfert 2s. 
Spectropolarimetric studies of Seyfert 2s (Miller & Antonucci 1983) showed the existence of broad 
emission lines also in some Seyfert 2s, which led to the formation of the unified model (UM) of 
AGN (Antonucci & Miller 1985). According to this model the same engine is at work in all kind 
of Seyfert galaxies, and differences between Seyfert Is and Seyfert 2s are due to the presence of 
an anisotropic absorber, which covers the broad line region (BLR) in Seyfert 2s. This absorber 
is often associated to a molecular torus. Several other lines of evidence support the UM, like the 
deficit of ionizing photons in Seyfert 2 galaxies, which indicates that the ionizing source is hidden 
from direct view (Schmitt & Kinney 1996), and the fact that the Narrow Line Regions (NLRs) 
are seen with conical shapes in Seyfert 2 galaxies and halo-like shapes in Seyfert Is (Pogge 1988). 
However, in the last years evidence of differences between Seyfert Is and Seyfert 2s has been dis- 
covered. Seyfert Is with significant absorption have been found (e.g., Cappi et al. 2006), along 
with Seyfert 2s without X-ray absorption (e.g., Bianchi et al. 2008), and Seyferts with significantly 
clumpy absorbers (e.g., Ricci et al. 2010). Furthermore spectropolarimetric surveys indicate that 
only ~ 30 — 50% of the Seyfert 2s show polarized broad lines (PBL), which might imply that not 
all of them harbor hidden BLRs (Tran 2001,2003). The hard X-rays are particularly well suited to 
test the unified model, as at these energies the photons are unaffected by photoelectric absorption 
and it is possible to have a direct view of the X-ray source, at least up to values of the column 
densities A^h — a? • For Ah > Gj 1 Compton scattering starts in fact to play an important role, and 
only < 50% of the incoming photons are detected. Thus, considering only Compton-thin objects, 
one would expect to detect consistent average spectra for different classes of AGN. 

2. Sample and data analysis 

The sample we used consists of all the 172 Seyfert galaxies detected by INTEGRAL IBIS/ISGRI 
(Ubertini et al. 2003) at z<0.2 during its first 7 years of observations. Of these, 45 are Seyfert Is, 31 
Seyfert 1.5s, 67 Compton-thin Seyfert 2s, 11 Narrow Line Seyfert Is (NLSls), 11 Compton-thick 
(CT) Seyfert 2s and 9 LINERs. The optical classifications were taken from Veron-Cetty & Veron 
(2010). To extract the average spectra of the different AGN samples we followed the procedure 
adopted by Walter & Cabral (2009). We created 500 x 500-pixels mosaic images modifying the 
coordinate system of each individual image, setting the coordinates of each source of the sample 
to an arbitrary fixed position (a=0, 5=0). The geometry of the image was also modified to have 
consistent PSF whatever the position of the source in the field of view (FOV). These mosaic images 
provide a stack of all the selected IBIS/ISGRI data for each considered sample. Individual sky im- 
ages for each pointing were produced in a broad energy band (i.e. 17-250 keV), and divided in 10 
bins. The spectra were extracted from the mosaics using mosaic_spec. In ngure[I]we show, as 
an example, a composite of the central part of the 17-80 keV mosaic image obtained for Seyfert Is. 
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3. Absorption in the X-rays 

Both photoelectric absorption and Compton scattering must be taken into account when mod- 
eling the effects of absorption in the X-rays. The photoelectric cross section o p h has a strong 
dependence on the energy, whereas Compton scattering depends on the Thomson cross section Oj, 
and it is constant with the energy, at least up to the Klein-Nishina decline. The ratio between Klein- 
Nishina and Thomson cross sections (7k n/<?t starts to diverge from the unity at E ~ 20 keV, and 
it is of ~ 0.6 at 200 keV. Compton processes become significant for column densities higher than 
(Tj 1 = 1.5 x 10 24 cm~ 2 . The cumulative effect of the two cross sections is given by: 

M(E) = e - a " h(E)Nn x e - a ? N *. (3.1) 

In Fig.|2| we show the effect of absorption (considering both photoelectric absorption and 



Compton scattering, as reported in Eq. 3. 1 ) on a power law model with a photon index of T = 1 .95 
in the 0. 1-300 keV energy range for different values of Ah- In Fig. || we show the fraction of 
escaping flux for different values of Ah in the energy band we used. About 60% of the original 
emission is unabsorbed for Ah = 7 x 10 23 cm" , and we used this value as a threshold between 
Compton-thin and Compton-thick sources. In our final sample we considered only AGN with a 
value of the column density Ah < V x 10 23 cm~ 2 . The values of the column densities were taken 
from soft X-ray (E < lOkeV) observations. We excluded from our sample the sources for which 
no value of Ah was available in the literature. 



4. Model-independent spectral analysis 

The hard X-ray spectrum of AGN can be well described by a power law continuum with an 
exponential cut-off, and a reflection hump peaking at E ~ 30 keV (Magdziarz & Zdziarski 1995). 
One of the drawbacks of a model dependent analysis is the parameter degeneracy. This might not 
allow to well constrain several parameters at the same time. An alternative method to characterize 
differences and similarities between different classes of Seyfert galaxies, independently of their 
average flux, is through a model independent approach. This has been done by normalizing the flux 
of the different spectra at the first bin (i.e. 17-22 keV), and then calculating their ratio. In Fig. || 



3 



INTEGRAL and the unified model of AGN 



Claudio Ricci 



10" 


i 1 1 


1 1 




10"' 




INTEGRAL IBIS/ISGRI band 




















N^ = 10 21 cm" 2 










N^ = 10 22 cm" 2 




'0 1 






N^ = 10 2i cm" 2 










N^ = 7x1 25 cm" 2 

N^ = 1.5x10" cm" 2 

IC=5x1Q 14 cm" 2 




1 0" 5 
10" 4 








10" 5 


r / / /' 








ID" 6 


= \ ; / // 

' i' ; / / 






in- 7 


1 ■ , i , / 1 


, , 1 ^ 





Energy (k.V) 

Figure 2: Effect of photoelectric absorption and Compton scattering on a power law model with a photon 
index of V = 1 .95 in the X-rays. 
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Figure 3: Effect of photoelectric absorption and Compton scattering on a power law spectrum (T = 1 .95) in 
the hard X-rays. The upper panel shows the observed spectrum for different values of the hydrogen column 
density A^h, and the lower panel shows the fraction of escaping flux in the 10 energy bins used in this 
work. In the lower panel squares represent the case N H = 10 23 cm- 2 , stars %=7x 10 23 cm , diamonds 
JV H = 1.5 x 10 24 cm- 2 and triangles iV H = 5 x 10 24 cm- 2 . 

we show the spectra and ratios obtained comparing the average hard X-ray spectrum of Seyfert Is 
to those of Seyfert 1.5s and NLSls. As it can be seen from the figure the average spectra of both 
Seyfert 1.5s and NLSls are consistent with that of Seyfert Is. NLSls are known to have steeper soft 
X-ray spectra than Seyfert 1 and Seyfert 1.5 galaxies. The fact that our sample of NLSls is hard 
X-ray selected probably introduces a bias, being the objects with a flatter spectrum more easily 
detected than those with a steeper spectrum. Malizia et al. (2008) studying a smaller sample of 
NLSls found that their hard X-ray spectrum is on average steeper than that of Seyfert Is. The 
difference between our results and theirs is probably due to the different sample, lower exposures 
and smaller energy ranges used in their work. 

Given the similarity of the spectra of Seyfert Is and Seyfert 1.5s, we re-ran the analysis consid- 
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ering all the objects of the two samples, obtaining an average spectrum of Seyfert Is and Seyfert 1.5s. 
We compared this spectrum to that of Compton-thin Seyfert 2s, obtaining a result remarkably dif- 
ferent from what we had for the other classes. As it can be seen from Fig. ||, the average hard 
X-ray spectrum of Compton-thin Seyfert 2s is significantly harder in the 22-60 keV band, showing 
a bump which might be linked to a stronger reflection hump. At energies higher than 60 keV the 
ratio of the two spectra is consistent with the unity, which implies that, although showing a stronger 
reflection component, Compton-thin Sy 2s have the same continuum of Sy Is and Sy 1.5s. 
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Figure 4: Spectra (upper panels) and ratios (lower panels) between the normalized spectra of Seyfert Is and 
Seyfert 1.5s (left), Seyfert Is and NLSls (right). 
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Figure 5: Ratio between the spectra of Compton-thin Seyfert 2s and the one of Seyfert Is and Seyfert 1.5s. 



5. The unified model of AGN 

The continuum emission of different classes of Seyfert galaxies has on average the same char- 
acteristics, which confirms the idea that different classes of objects have the same engine, as pre- 
dicted by the zero-th order unified model of AGN. The stronger reflection observed in the spectrum 
of Compton-thin Seyfert 2s cannot be explained in terms of absorption, as the chosen threshold of 
A^h should guarantee a transmission dominated spectrum. Moreover considering that the average 
column density of our sample of Seyfert 2s is Nu = 1 .5 x 10 23 cm~ 2 , one would expect only ~ 10% 
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of the flux to be Compton scattered. The excess in the 20-60 keV band might be related either to 
an intrinsic stronger reflection component of Seyfert 2s, or to the presence of clumps in the torus, 
as predicted by Elitzur & Shlosman (2006). 

6. Conclusions 

We presented the study of the average INTEGRAL IBIS/ISGRI hard X-ray spectra of different 
classes of Seyfert galaxies. We found that, as predicted by the UM, the average spectra of Seyfert Is 
and Seyfert 1.5s are consistent in the 17-250 keV band. NLSls also show a spectrum similar to 
that of Seyfert Is and Seyfert 1.5s, probably due to the hard X-ray selected nature of our sample. 
The average emission of Compton-thin Sy 2s is harder than that of Sy Is and Sy 1.5s in the 20-60 
keV energy band, which might be due either to an intrinsically stronger Compton reflection or to 
a clumpy torus. The UM is confirmed on the zero-th order, although an homogenous anisotropic 
absorber alone cannot easily explain the differences between the spectrum of Seyfert 2s and that of 
Seyfert Is and Seyfert 1.5s. 
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